Another Origin of Bimodal Size Distribution in InAs Self- Assembled Quantum Dots 
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The evolution of InAs quantum dots grown on InP substrates by metal-organic vapour phase 
epitaxy is studied as a function of InAs coverage. Under specific growth conditions, the onset of 
the two- to three-dimensional transition is seen to proceed via two distinct pathways: through (i) 
an abrupt appearance of quantum dots as expected in the usual Stranski-Krastanov growth picture 
and (ii) a continuous evolution of small surface features into well developed quantum dots. The 
average size of the features in both these families increases with coverage, leading to a bimodal 
distribution in dot sizes at an intermediate stage of growth evolution that eventually becomes a 
unimodal distribution as more material is deposited. Complementary information obtained from 
independent measurements of photoluminescence spectra and surface morphology is correlated and 
is found to be independently consistent with the picture of growth proposed. 

PACS numbers: 68.65.Hb, 78.67.Hc, 81.07.TVi, 68.37.Ps, 81.15.Gh 



The semiconductor heterostructures grown via the 
strain-mediated Stranski-Krastanov (S-K) routeQ, Q 
have been of special interest. The resulting self assem- 
bled clusters are extremely small (~10nm, hence called 
quantum dots), coherently strained (hence optically ac- 
tive) and have a size dispersion that may be acceptable 
for their use in real optoelectronic devices, e.g., 0, 
The S-K transition has been systematically studied in 
most III-V semiconductors of interest (e.g. InAs/GaAs 
H IS II lljbi-sGaxAs/GaAslIft E5, InAs/InP[lp, 
GaSb/GaAsfl3j|). While the growth route may be ratio- 
nalized by energetic considerations that predict an ini- 
tial two-dimensional growth followed by an abrupt ap- 
pearance of three-dimensional clusters upon the 'wetting 
layer'[l|, the actual process of self-assembly is compli- 
cated by kinetic effects, substrate conditions and ori- 
entation, a non-quiescent wetting layer during growth 
evolution, alloying and long-ranged substrate-mediated 
elastic interactions. For example, it has also been ob- 
served, though not consistently by different groups, that 
there may be redistribution of matter leading to a de- 
crease in the quantum dot (QD) size with coverage. Of- 
ten a bimodal size distribution |S ^3 has also been ob- 
served at an intermediate growth stage. Although many 
of these phenomena have been described in terms of elas- 
tic energy barriers corresponding to change in shape and 
a maximum permissible size for quantum dots beyond 
which they are dislocated, the complex growth evolution 
of self-assembled dots, especially under far from equilib- 
rium (usual) growth conditions, is still very much an open 
problem despite more than a decade of extensive work. 

In this paper, we focus on InAs/InP self-assembled 
QDs grown by metal-organic vapour phase epitaxy 
(MOVPE)0 G]|. The samples are studied at differ- 
ent stages of evolution by independently studying and 
quantitatively correlating the morphology and the optical 
emission spectra. Contrary to many of previous reports, 
we have, on samples grown under more non-equilibrium 
conditions (lower temperature, 450° C and higher growth 



rate ~l-2 monolayers (ML) per second) where surface ki- 
netics may be an important limiting step, observed that 
the 2D-3D transition is not completely abrupt and pro- 
ceeds via two different pathways. This study also yields 
a different picture for bimodality in the QD's size distri- 
bution. 

Growth was carried out using low pressure (100 torr) 
MOVPE on n+ doped (001) "epi-ready" InP substrates 
in a horizontal reactor with hydrogen as the carrier gas. 
Group III and V sources were trimethyl-indium and ar- 
sine and phosphine respectively. Prior to InAs growth, a 
InP buffer layer was grown first at 625°C and then with 
the temperature continuously ramped down and stabi- 
lized to 450° C. InAs layers were grown at a relatively 
low temperature of 450° C at a growth rate0| of ap- 
proximately 1.8ML/s. A pair of samples was grown with 
identically deposited InAs layer in two growth runs. In 
the first case, the sample was taken out of the reactor 
after InAs deposition itself to enable a study of surface 
morphology and in the second case a InP cap layer was 
grown for samples used for photoluminescence study. For 
these samples, about 50 A InP was deposited at the 
InAs deposition temperature to avoid any further ripen- 
ing during the higher temperature overgrowth. Surface 
features were studied ex-situ using an atomic force mi- 
croscope (AFM) in contact mode, typically within a few 
hours of sample growth. The PL spectra were measured 
at ~25 K at low enough excitation power (~ 0.5W/cm 2 ) 
to preclude any subband filling and were corrected for 
the system response using a standard black body source. 

Fig. U (inset) shows \[im x 1/im AFM scans for sam- 
ples with progressively increasing InAs coverage in Fig. 
^a) to (d). The histograms of the island heights corre- 
sponding to these AFM images are shown in the main 
Figffl All well-separated (^20nm) convex features with 
heights above two ML were counted here. Fig^ (a) de- 
picts an early stage of growth. Here we observe that 
along with very few (~ 5 x 10 8 cm' 2 ) well developed 
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FIG. 1: Histogram of island heights inferred from 1/im x ljitm 
AFM images on samples with increasing InAs coverage. The 
coverage is approximately: (a) 4-5ML, (b) 6.5ML (c)9ML 
(d)14-5ML. The coverage is calibrated against the WL PL peak 
in FialWa). See text. The dotted line depicts the separation 
of the islands into two families, evidently evolved from two 
different routes (see Fig. 3). This leads to a bimodal size dis- 
tribution (a)-(c) which merges into a single broad distribution 
(d). 



(height > 8 run) QDs, there are many small surface fea- 
tures of minimum height of ~ 14A and a mean height 
of 26A contributing to a large number of counts at low 
heights. Furthermore, in Fig. |U(b)-(c), we observe that 
these small surface features are stable and continuously 
grow in size as more material is deposited. Therefore, 
with incremental coverage, Fig^b)-(d), the morphology 
evolves as (i) matter is accreted by all pre-existing fea- 
tures making them grow in size (ii) new well developed 
QDs of height ~ 8nm are spontaneously generated (usual 
S-K transition). Growth via two independent pathways 
naturally leads to a bimodal size distribution at an inter- 
mediate stage of growth, Fig^c). The role of existing 
(two- and quasi-three-dimensional) surface features dur- 
ing the S-K transition has been much discussed [(| 0,0]. 
Here these surface corrugations are only being connected 
to dots which continuously evolve in size and we have 
not examined their role, if any, as precursors to abruptly 
generated well developed dots (S-K transition). 

The above mentioned picture of growth evolution is in- 
dependently observed in optical properties. FigEI shows 



FIG. 2: (Solid line) 25K PL Spectra. The growth conditions 
for samples whose response is depicted in (a)-(d) are the same 
as their counterparts in Fig. QJ except for the InP overgrowth,. 
(Squares) Calculated transition energies scaled with the 
heights histogram (Fig. QJ) in an attempt to reconstruct the 
PL spectra. The expected response from only the features 
higher than the cut-off depicted by dotted lines in Fig. Q is 
shown. Note that there is no fitting parameter. The higher 
energy feature is fitted to two Gaussians. (Dotted line) Wet- 
ting layer contribution to PL. PL from the wetting layer ap- 
pears at 1.02—1.03 eV in (a)-(c) and shifts to 1.10 eV for 
(d). See text. The other fit corresponding to the evolution of 
small surface features is not shown here for brevity. 



the low temperature PL spectra measured on the four 
similarly grown but capped samples. Along with the 
characteristic PL from well-developed dots, the PL spec- 
trum also shows a wetting layer peak around 1.02—1.03 
eV that is broadened towards lower energy by the emis- 
sion from these optically active quasi-three-dimensional 
clusters. The spectrum in Fig[5{c) shows almost three 
distinct peaks corresponding to the wetting layer and the 
two kinds of quantum dots. Finally in the saturation 
regime (when the number of well resolved QDs does not 
increase with coverage), the distribution of the dot sizes 
is again roughly unimodal, but with a large dispersion. 

The optical signature from the wetting layer was in- 
ferred by fitting two Gaussians, (corresponding to signals 
from WL and quasi-three dimensional clusters) to the 
high energy PL peak. We observe that the wetting layer 
peak is constant to within lOmeV in Fi^fa)-(c). This 
low temperature PL peak at around 1.02— 1.03eV corre- 
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sponds to an approximately 4.5 ML strained InAs/InP 
quantum well|l8j. This value may be used to calibrate 
the thickness with the deposition time, which is other- 
wise difficult to estimate with sub-ML precision in a typ- 
ical MOVPE set-up without in-situ diagnostic tools. It 
is worth pointing out that there has been a wide vari- 
ation in the reported value of the wetting layer thick- 
ness for InAs/InP, with the lower limit being ^1ML 
0,H(j. With a lower growth temperature, the wetting 
layer thickness is expected to be (exponentially) larger 
than its equilibrium value |2l| . since the former is in- 
versely dependent on the diffusion length [22|. Because 
450° C is among the lower reported growth temperatures 
for InAs dots, a thicker wetting layer is naturally ex- 
pected. Remarkably, a recent reflection high energy elec- 
tron diffraction study p3| on 2D-3D transition in MBE 
grown InAs/InP films, reports for growth at 450°C, a 
value of the wetting layer thickness that is very similar 
to what we have inferred from the PL spectra. 

Finally, in Fig. Etd), the wetting layer PL shifts to a 
higher energy, ^1.10 eV indicating the well known nar- 
rowing of wetting layer in the saturation regime [24| . 

The observed PL spectra from QDs should in prin- 
ciple be derivable from the quantitative analysis of the 
AFM images |25l. Using published ground state energy 
calculations [2^for InAs/InP QDs, the expected emission 
spectra from large dots are shown in FigEJopen squares). 
While there is a good quantitative agreement without any 
fitting parameter for the expected and observed spec- 
tral features from large dots (those above the size cut-off 
depicted by dotted lines in Fig. QJ, the PL from the 
smaller features in each histogram did not agree as well 
and is therefore not shown in FigEl It has been previ- 
ously established that the capping process can consider- 
ably change the morphology of the individual quantum 
dots and it is likely that the smaller dots are more af- 
fected by capping [2?], HI] and compositional changes due 
to the As/P exchange during the early stage of InP over- 
growth. 

Implications: We have observed that the 2D-3D transi- 
tion is not abrupt and the bimodality in the QD's size 
distribution at an intermediate stage of growth is a conse- 
quence of the growth evolving via two distinct pathways. 
These two observations are contrary to previous studies 
on InAs/GaAs @, Q and In^Gai-^As/GaAs ^(| where 
the 2D-3D transition had been found to be abrupt and 
thermodynamically first-order with the surface coverage 
playing the role of a critical parameter 5J. An abrupt ap- 
pearance of quantum dots beyond a critical coverage has 
also been theoretically reproduced within a rate equation 
based model |2^| . 

Bimodality in dot sizes also implies that the QDs en- 
semble cannot be characterized by a single length scale 
(e.g., mean island size)Q during most of the evolution 
and rules out the very attractive possibility of data col- 
lapse onto a universal scaling.3C} function during inter- 
mediate stages of growth. Absence of scaling at an in- 
termediate growth stage has actually been observed by 



Krzyzewski, et al. (e.g., see ref. p}). The issue of bi- 
modality in dot sizes at intermediate coverage, although 
not fully resolved, is typically explained by the presence 
of energy barriers |3l) corresponding to another (shape- 
change, from pyramid to dome) first order transition 
that has been experimentally studied in SiGe/Si [32| . 
GaN/AIN |H and InAs/GaAs 
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FIG. 3: Variation of the island density with coverage for the 
two families of islands corresponding to the two modes in the 
heights distribution in Fig. 1. Solid square represents the late 
stage when the distribution is umimodal (Fig. YWd) ). 

It is instructive to (visually) separate § the heights 
histograms in Fig. ^ at the minima between the two 
modes (dotted lines in Fig|TJ. Figgis the corresponding 
plot of cluster densities in these two families (S-K dots 
and continuously evolving clusters) as a function of cover- 
age. The increase in the density of larger dots with cover- 
age is abrupt, whereas the number of smaller features de- 
creases with coverage due to coalescence. At large enough 
coverage the distribution becomes unimodal. Without in- 
voking the existence of energy barriers, this is most sim- 
ply understood as being due to the difference in growth 
rates of larger and smaller dots (due to their different 
volumes). 

Among the various studies of growth evolution, the 
results in references 0, |J, although on MBE grown 
InAs/GaAs, are qualitatively most similar to ours. Nev- 
ertheless there are some important differences both in 
data and interpretation of results. Firstly, the quasi- 
three dimensional clusters were thought to be precursors 
to all the quantum dots giving a common origin to the 
dots in both the families. Possibly due to a smaller strain 
in the InAs/InP system as compared to InAs/GaAs, we 
observe an uninhibited increase in size and a correspond- 
ing redshift in PL with coverage. This is in contrast 
with InAs/GaAs QDs, where a pronounced barrier seems 
to restrict the maximum dot size to ~ 8nm. This sup- 
ports the general observation of larger dispersion in dot 
sizes and low temperature PL linewidths (~ lOOmeV) in 
InAs/InP. 
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The differences observed in growth evolution by differ- 
ent groups may be a result of differing substrate/buffer 
layer conditions. It is possible for the local roughness in 
the wetting layer surface to stabilize pre-existing quasi- 
three dimensional surface structures and determine the 
extent of bimodality at an intermediate growth stage. 
Many of the thermodynamic arguments used to describe 
the QD self-assembly are relevant only under quasi- 
equilibrium conditions of growth which many such pre- 
vious experiments tried to maintain. Our results, on 
the other hand, are obtained under much higher growth 
rates, more typical growth conditions of quantum dots. 
Conclusions: Studying the growth evolution of MOVPE 



grown InAs/InP self-assembled quantum dots, we have 
observed that an alternate pathway for the 2D-3D tran- 
sition exists that naturally explains the often observed 
bimodality in the quantum dots size distribution at an 
intermediate growth stage. Independent evaluations of 
the morphology and the photoluminesence spectra are 
consistent with the picture of growth presented. 
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